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Single cell mechanics: stress stiffening and kinematic hardening 
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Cell mechanical properties are fundamental to the organism but remain poorly understood. We 
report a comprehensive phenomenological framework for the nonlinear rheology of single fibroblast 
cells: a superposition of elastic stiffening and viscoplastic kinematic hardening. Our results show, 
that in spite of cell complexity its mechanical properties can be cast into simple, well-defined rules, 
which provide mechanical cell strength and robustness via control of crosslink slippage. 

PACS numbers: 87.15.La, 83.60.Df, 83.60.La, 87.16. Ka 



Intracellular transport, cell locomotion, resistance to 
external mechanical stress and other vital biomcchani- 
cal functions of eukaryotic cells are governed by the cy- 
toskeleton, an active biopolymer gel [lj . This gel is made 
of three types of biopolymers, actin, microtubules and 
intermediate filaments, crosslinked by a multitude of pro- 
teins with different properties in terms of connection an- 
gles, bond strengths and bond lifetimes. The actin cy- 
toskcleton -the major force-sustaining structure in our 
experiments- is made of filaments of about micrometer 
length and presents a weak local structural order. The 
cytoskeleton also comprises molecular motors, proteins 
that move on actin or microtubule filaments driven by 
chemical energy. How the cytoskeleton in conjunction 
with biochemical regulatory circuits performs specific, 
active mechanical tasks is not understood. When cells 
attach to biological material they often biochemically rec- 
ognize the binding partner. The cytoskeleton organizes 
accordingly and produces a mechanical response. Ac- 
tive cell responses such as contraction are well separated 
from passive rheological properties by their timescales Q • 
Passive rheological cell properties have been studied with 
various techniques on subcellular and supcrcellular scale 
Q . From the measurements with different techniques on 
different eukaryotic cell types a broad relaxation spec- 
trum arises as a common feature of passive, linear cell 
rheology 0, H[ . The description of the non-linear regime 
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In the following we present microplate rheology exper- 
iments where individual cells are stretched between two 
plates (Fig. [lj . The advantage of the setup is that the 
cells possess a well-controlled geometry and adhere via 
chosen biochemical linkers, which better define the cy- 
toskclctal state. Quasi-differential cell deformations re- 
veal an elastic stiffening response. The corresponding 
nonlinear clastic modulus depends on the cell prestress 
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Figure 1: A fibroblast adhering between two microplates. a: 
Right after contact. Bar: 10 /im. b: After ~20 min at 35°C, 
strongly adhering cells often adopt a concave shape. From 
the apparent diameter Do (arrows) we estimate the initial 
area Aq := n(Do/2) 2 . c: Under large stretch. 

but is independent of cell length. Large deformations re- 
veal an inelastic regime with a (counterintuitive) linear 
force—length relation. Both relations simply superpose 
to generate the response to more complex deformations. 
The cell response reduces to the integral of the differential 
measurement when the inelastic response is abolished by 
fixation. Hence, in spite of the complexity of the eukary- 
otic cell cytoskeleton and large cell heterogeneity, passive 
nonlinear cell rheology can be reduced to simple rules. 

Experimental setup. — We refer to 0, 0, EH f° r 
details. A 3T3 fibroblast [l3| adheres between two 
fibronectin-coated glass microplates. One of them is flex- 
ible: its deflection gives the perpendicular force F acting 
on the cell. This plate is translated by a piezoelectric 
actuator controlled by a personal computer. The com- 
puter calculates force F and cell length £, and adjusts 
the piezo position via a proportional feedback loop to 
impose a given experimental protocol. Experiments are 
performed at 35°C, in standard medium with Lysophos- 
phatidic acid 50 fjM (Sigma). Cells are left to adhere 
for 30 min before measuring. All results described here 
are fully reproducible for fibroblasts which adhere suffi- 
ciently strongly to sustain pulling forces of 100 nN for 
several hours, which means about 30% of the cells in cul- 
ture. 

Loading and unloading at constant rate — Wc 

stretch the cell by 100% at a constant rate I while mea- 
suring the force F (Fig.[2Js.). The slope dF/d£ initially de- 
creases, reaching a constant value at an elongation ~ 10% 
(Fig. [lb). Beyond 10% and up to 100% stretch, the F(£) 
relation is in most cases linear. After loading, the cell 
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Figure 2: Loading and unloading at constant rate £. a: length 
I and lagrangian stress F/Ao (where Ao is the initial area) as 
a function of time. After each ramp, F stabilizes at a nonzero 
value 5 (dotted line), b: F/Aq as a function of £/£o (where £o 
is the steady zero-force length) during loading and unloading. 
Black curve: £ = 1 /im/s. Grey curve: i — 0.03 /im/s. 




length I is held constant for a few minutes; the force F 
relaxes to a steady non-zero value J which does not evolve 
faster than ~ 1 nN/s. An analogous response is observed 
upon unloading. The procedure is repeated with differ- 
ent rates I between 3 nm/s and 10 /im/s. The asymptotic 
slope dF/d£ and the equilibrium force 3 arc independent 
of the loading rate in the explored range. 

Small amplitude stiffening, large amplitude 
linearity — To explore small and large deformation 
amplitudes simultaneously we perform a loading ramp 
with superimposed harmonic oscillations, imposing 
£(t) = vt + sin(w<). Fig. [3] shows a typical ex- 
periment. The response to small oscillations indeed 
stiffens with increasing stress. Yet, the averages over 
an oscillation period of the force (F) and length 
{£) are linearly related as inferred from the position 
of the loops. Therefore, we observe both responses 
simultaneously: stiffening at small amplitudes, as 
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Figure 3: Small amplitude stiffening, large amplitude linear- 
ity, a: F and £ as a function of time, b: F/Ao as a function 
of £/£o. For clarity, only a few loops are shown. The dashed 
line highlights the linear relation between the average values. 
Inset: Differential modulus j0| as a function of average force 
(F) for the data shown in b. The response to the small oscilla- 
tions shows stiffening, following the master-relation reported 
in 0. c: Eqs. 1-4 taking 7 = 5 and V = 0.01 dF/d\. 



Figure 4: Elastic / inelastic behavior, a: Imposed length £ as 
a function of time, b: F/Ao as a function of £/£o for a given 
cell. Reversible (elastic) behavior upon direction reversal is 
observed only close to a previous turning point, as in C, E, H. 
The response becomes irreversible (inelastic) after a steady 
large deformation: at the turning points D, F, G, I, the F(£) 
curve does not retrace its previous path. Between F and I 
the experiment is equivalent to C-F, but in the unloading 
direction. The response is seen to be equivalent, showing 
the sense of deformation to be irrelevant, cl: Another cell, 
probed at a rate £ = 0.1 /im/s. c2: Same cell as cl but at 
£ — 1 /im/s. dl: Another cell, 0.1 /im/s. d.2: Same cell as dl, 
1 /im/s. e: Prediction of the constitutive relation (Eqs. 1-4). 

reported in Q, and linearity at large amplitudes. 

Small amplitude reversibility, large amplitude 
irreversibility — We study the amplitude dependence 
at a constant deformation rate \£\. An essential feature of 
the protocol (Fig. \4&) are the turning points separated at 
various distances in order to study the reversibility of the 
response. Similar procedures can be found in plasticity 
textbooks [lj]. As Fig.rjJ) shows, the reversibility of the 
response upon a change of direction is determined by the 
distance to the previous turning point. Where turning 
points are separated by less than 10% stretch, the re- 
sponse is reversible (elastic). More than 10% stretch be- 
yond a turning point, the response becomes irreversible 
(inelastic): the F(£) curve does not retrace its path upon 
direction reversal. In this inelastic regime the F(£) rela- 
tion is approximately linear. Its nonzero slope leads to a 
translation of the elastic region by the inelastic deforma- 
tion, a behavior known in plasticity as linear kinematic 
(or directional) hardening [3, EM EE E3 • Alternatively, 
the inelastic contraction under pulling tension between X 
and G in Fig. 0b is a strong Bauschinger effect (a decrease 
in yield stress upon unloading) [3, EH . 

Large amplitude stiffening after glutaraldehyde 
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Figure 5: a: Imposed length I as a function of time. A sin- 
gle loading cycle with amplitude 30% is performed, b: Force 
F as a function of length £. "Normal": before adding glu- 
taraldehyde. "Fixed": in presence of glutaraldehyde 0.1%. 
The dotted line is a fit to Eq[5] c: The derivative dF/d£ of 
the "fixed" curve (black line) plotted against the differential 
master-relation (gray dots, as discussed in Ref. (?]])• 



fixation. — We add glutaraldehyde 0.1% in order to pre- 
vent slippage of cytoskeletal connections. Loading at con- 
stant rate (Fig. [5ji) reveals a positive curvature d 2 F / d£ 2 
(Fig. [5}d). The numerical derivative of the F(£) relation 
obtained from fixed (hence dead) cells is the same as the 
differential master-relation obtained on living fibroblasts 
(Fig. [3b, from Rcf. >]). The F(£) relation after fixation 
closely resembles the exponential stress-stretch relations 
known from whole tissues [l9|, [2(j ■ 

Rate dependence— In the inelastic regime the width 
of the hysteresis loops increases with stretch rate (Fig. 
[6^i). To characterize this rate-dependence, we define 
the overstress AF as the extent of force relaxation af- 
ter unloading (Fig. [6^,). The overstress AF behaves as 
log(d£/dt), approaching zero at a non-zero pulling speed 
of about 10 nm/s (Fig|6fc>). Below such rates the behavior 
becomes active and erratic and the overstress ill-defined. 

Viscoplasticity— We now propose a minimal consti- 
tutive relation for fibroblasts under uniaxial extension. 
First we decompose the measurable cell length £ into 
inelastic rest length £ and elastic stretch ratio A, 

£ = XZ. (1) 

The force is a function solely of the elastic strain, 

F oc (A - i) e ^ 2 +2/A-3) : (2) 

where for concreteness we use exponential elasticity fl9| , 
according to Fig.[5fc>, c. As a flow rule relating the inelas- 
tic strain rate £ to the force F, we propose an exponential 
function of the overstress F — J, according to Fig. [5b>: 



£ = ?;sgn( J F-30el F ~ s l/ x ' 



(3) 



with a vanishing dissipation as the flow rate £ approaches 
v. The equilibrium force 5 is the essential internal vari- 
able to describe kinematic hardening 0, [lH, QjJ ■ The 
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Figure 6: Loading and unloading for several rates £. a: F as a 
function of £/£q for a given cell, where £q is the initial length. 
Rates are: 1) 10 nm/s, 2) 30 nm/s, 3) 0.1 /im/s, 4) 1 /im/s. 
The overstress AF is denned as the extent of relaxation after 
unloading, b: AF as a function of rate £ for 7 different cells. 



drag force T> sets the scale where the overstress induces 
inelastic flow. To obtain an increased dissipation at large 
forces (e.g. the increase in the area of the loops in Fig. 
[3]), we take the drag T> as proportional to the nonlin- 
ear modulus, dF/dX. Finally, we have linear kinematic 
hardening: 
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(4) 



This is an empiric description along the lines of modern 
viscoplasticity [IJ, [l?], [Hj], without explicit history de- 
pendencies. As Figs. [3b and [4^ show, it captures the 
essence of the phenomenology. At small amplitudes, 
\F — $\ -C T>, the deformation is essentially elastic: I ~ A. 
At large amplitudes the overstress \F— $\ approaches the 
drag stress V and the deformation becomes increasingly 
inelastic: £ ~ £ oc $ ~ F. Nevertheless, this constitutive 
relation is not yet a full description. The details of the 
linear regime [j, [ll| and fluidization at large flow rates 
[?J still have to be incorporated to it, whereas active con- 
traction and inelastic deformation at rates below v, as 
well as the force fluctuations seen in Fig. [2] may require 
a different approach. 

Discussion — Our glutaraldehyde fixation experi- 
ments show that stress stiffening in fibroblasts Q is due 
to the nonlinear elasticity of the cytoskeleton, unrelated 
to biological signalling or restructuring. In agreement, 
very similar stiffening is known from biopolymer net- 
works [2l|. To date the precise microscopic mechanism 
remains unclear; stretching [U and bending % [23| 
of single filaments as well as filament alignment [2J] have 
been proposed. 

The energy to reach the linear inelastic regime in 
loading experiments (from Fig. [3]) is — - IfiN (im ~ 
2.5 10 8 ksT. During stretch, first elastic elements must 
be loaded, until they dissipate the stored energy upon 
bond rupture. Taking typical orders of magnitude for 
bond energies [25| and dissipation [26| and a mesh size 
of 100 nm [l|, over 10% of the actin cytoskeletal bonds 
must be ruptured to reach stationary flow. However, in 
order to observe stress stiffening ubiquitously during the 
inelastic deformation, the actin gel must be always above 
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percolation threshold. For a gel with one bond relaxation 
time one expects a rate-independent overstress [27j |: we 
observe an exponential ratc-overstress dependence. In 
general, dissociation rates depend on the for ce p er bond 
/as - e f / fo , with a force scale / ~10 pN [1| 11]. In 
agreement, in fibroblasts the drag T> is about 100 nN, cor- 
responding to 1-10 pN per filament. Interestingly, the 
inelastic stretch rate where the cell flows without hys- 
teresis is of the order of 10 _3 s _1 , a typical rate for active 
processes such as crawling and contraction [l|, Q . Thus, 
spontaneous bond dissociation may be what limits active 
phenomena to long timescales 

A living cell can neither be purely elastic, nor possess 
a yield stress within the physiological "working range" . 
Kinematic hardening viscoplasticity can be understood 
as a consequence of these conditions. Combined with 
a sharply rising rate-ovcrstress dependence, it prevents 
cell breakage in our large deformation experiments: a 
cell portion under increased stress will readily flow and 
increase its equilibrium stress to reach a stable situation. 
Rather than break at a given spot, the cell prolongs ho- 
mogeneously along its length. This homogeneous defor- 
mation may be behind the robust linearity of the kine- 
matic hardening response, since integration of a constant 
magnitude along the cell length naturally gives a linear 
length-dependence. However, identification of the precise 
molecular mechanisms behind this unusual behavior in a 



soft system is a task for the future. At least, one can rule 
out a trivial interpretation in terms of a Hookean spring 
element (in form of intermediate filaments, for example) 
in parallel with a stiffening viscoelastic liquid: since the 
liquid cannot sustain an average stress, a non-mechanical 
coupling between the two elements is needed for the cell 
to stress-stiffen. Thus, in a mechanical interpretation 
hardening and stiffening must originate in one and the 
same mechanical element. If intermediate filaments [l[ 
play a role, they must be interconnected with actin into 
a single network. This scenario reminds of composite 
alloys [HI and granular materials [l6|, where kinematic 
hardening arises as the inelastic flow induces alterations 
of directional nature to the microstructure. 

Summarizing, we have shown that cell mechanical 
properties in uniaxial stretching experiments can be thor- 
oughly described by the superposition of two simple re- 
lations: exponential elasticity, and viscoplasticity with 
linear kinematic hardening. Given the cytoskeletal com- 
plexity, this is unexpected. A complete picture of passive 
cell rheology spanning from molecular details to a simple 
phenomenological description and straightforward theo- 
retical concepts seems in close reach. 
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K. Kroy for inspiring discussions and support. This work 
has been funded by the Universitat Bayreuth. 
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